Spin Splitting and Spin Current in Strained Bulk Semiconductors 
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We present a theory for two recent experiments in bulk strained semiconductors P, 0| and show 
that a new, previously overlooked, strain spin-orbit coupling term may play a fundamental role. We 
propose simple experiments that could clarify the origin of strain-induced spin-orbit coupling terms 
in inversion asymmetric semiconductors. We predict that a uniform magnetization parallel to the 
electric field will be induced in the samples studied in P, |3 for specific directions of the applied 
electric field. We also propose special geometries to detect spin currents in strained semiconductors. 

PACS numbers: 72.25.-b, 72.10.-d, 72.15. Gd 



Spin manipulation in semiconductors has seen remark- 
able theoretical and experimental interest in recent years 
with the advent of spin-electronics and with the realiza- 
tion that strong spin-orbit coupling in certain materials 
can influence the transport of carriers in so-called spin- 
tronics devices 3]. In particular, the issue of creating 
spin polarization of carriers in nonmagnetic semiconduc- 
tors with spin-orbit coupling using only electric fields has 
caused a flurry of theoretical and exp erimental activity 

miBSBIiESIillilElElElEl two kinds of 

theories of spin-polarization under the action of an elec- 
tric fleld have been put forward. The first kind, dating 
back since the mid 1980's predicts the existence of a 
spatially homogeneous net spin polarization perpendic- 
ular to the applied electric current in two dimensional 
samples with spin-orbit interaction. This effect is dis- 
sipative and has been recently observed experimentally 
[itI . There also exist two very recent 0, |^ theories pre- 
dicting non-dissipative, intrinsic spin currents with the 
spin polarization and flow direction perpendicular to each 
other and to the electric fleld. This effect does not cre- 
ate a bulk magnetization but, if observed, can be used 
for spin injection, and its validity is being experimentally 
tested at the present time. One of the theories pre- 
dicts a spin current polarized out of plane and flowing 
perpendicular to the in-plane electric field applied on a 
2-dimensional semiconductor sample exhibiting Rashba 
spin-orbit coupling. As long as the Rashba spin split- 
ting is large enough, the spin conductivity is 'universal' 
{e/Swh) in the sense that it does not depend on the value 
of the coupling. The other effect P| appears in the va- 
lence band of the bulk samples and is proportional to the 
spin-orbit splitting of the valence bands (to the difference 
between the Fermi momenta of the heavy and light-hole 
bands). 

In the first part of this letter we analyze the theory be- 
hind two recent experiments in bulk strained semiconduc- 
tors P, 0] where an electric-field-induced uniform homo- 
geneous spin polarization upon an applied electric field 
is observed. We make the case that the observed spin- 
splitting (whose origin is puzzling) and spin polariza- 
tion is due to a previously overlooked strain-spin-splitting 



term, and propose easy experimental checks of our the- 
ory. 

In the second part of this letter we predict the ap- 
pearance of an intrinsic spin polarized spin current in 
n-doped bulk (and 2 dimensional) strained semiconduc- 
tors (GaAs, GaSb, InSb, InGaAs, AlGaAS, etc) under 
the influence of an electric field. The spin conductance 
is 'universal', in the sense that it does not depend on the 
value of strain (for large enough strain), but it is propor- 
tional to the average Fermi momentum of the conduc- 
tion band. The effect is due to the spin-orbit splitting of 
the conduction band under strain and is hence absent in 
strain-free semiconductors. The very long spin relaxation 
time in the conduction band as well as the relative pen- 
etration of strain engineering in semiconductor industry 
applications make this effect of potential technological 
importance. We propose an experimental technique us- 
ing the already existing setup in 0, |3| to measure the 
spin current and to differentiate between the intrinsic 
spin current and the uniform magnetization effects. 

In Pi nine samples of n-doped (n = 3 x 10^^cm~^) 
lUxGai-xAs {x = 5% — 7%) of thicknesses between 
200nTO — ISOOnm, grown in the [001] direction on un- 
doped GaAs substrate, are used to probe the electron 
spin dynamics through time and spatially resolved Fara- 
day rotation (FR). The length and width of the samples 
are roughly 300/im x 80/im. The lattice mismatch pro- 
vides for diagonal strain in the x,y,z = [100], [010], [001] 
directions of 0.04% — 0.46% 18] (contrary to claims in 
pSli] . the lattice constants in x and y directions are also 
strained, this being a generic feature of [001] growth). 
Moreover, anisotropic shear strain develops in all di- 
rections {xy,xz,yz = [110], [101], [Oil]) due to different 
direction-dependent strain relaxation rates at the growth 
temperature of around 500C This guarantees that 
all the components of the strain tensor e^, i,j = x,y,z 
are non-zero and of the same order of magnitude. The 
magnitude of the strain components is given in Table[l]. 

Pump-probe FR beams measure the total magneti- 
zation of the optically injected electron spins in the 
growth direction z when the samples are placed in an 
electric field on the [110] and [110] directions, respec- 



2 



tively. The dynamics of the spin packet is mainly de- 
scribed by a precession around a total magnetic field 
Btot — Bint + Bf,j.t where B^^t is an externally applied 
magnetic field whereas B^nt is the momentum-dependent 
internal magnetic field caused by the spin-orbit coupling. 
The precession around a Bint is the main feature of 
most of the spintronics devices, starting with the Das- 
Datta spin-field transistor . Under an applied electric 
field, the average particle momentum acquires a non-zero 
value, parallel to the electric field. The internal magnetic 
field is caused by the spin orbit coupling: the electric field 
acts on the particle momentum which in turn couples to 
the spin. The signal at the probe beam can be fitted to 
cos{g^B\Bint + Bextl^t/fi-) whcre is the Bohr magne- 
ton, g is the electron gi-factor while At is the temporal 
delay between the pump and probe pulses. This fit gives 
the direction and value of Bint which turns out to be per- 
pendicular to the applied electric field E and the z axis 
(for E in-plane); the value of Bint is used to determine 
the spin splitting Aq = g^sBint and a phenomenological 
relation Aq = l3vd is observed where Vd is the spin-drift 
velocity and /3 is a constant of proportionality that is the 
focus of the experiment 0|. Experiments find that Vd is 
linearly proportional to the electric field E. 

As a first step, let us theoretically address the ques- 
tion of origin of /3. By group theory, inversion symme- 
try breaking bulk strained semiconductors exhibit three 
main types of spin splitting ■ 
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Hi = X[aa:kx{kl - kl) + ayky{kl - kl) + (j,k,{kl - /c^)] 

H2 = '^C^\(7x{^xyky — ixzkz) + CFy{eyzkz — iyxkx)'^ 

~^ f^zi^zxkx ~ ^zyky)] 
H3 = D\axkx{^zz~^yy)~\''^yky{ixx^^zz)''rOzkz(,^yy~^xxY 



(1) 



TO = 0.0665too is the effective electron mass in the con- 
duction band ji^], X,C3,D > are material constants, 
^x.y^z are the 3 spin-Pauli matrices, and e^j, i, j — x,y, z 
are the components of the symmetric strain tensor. 

All three Hamiltonians can be written as the coupling 
of a fictitious fc-dependent internal magnetic field Bint 
to the electron spin, Bint{k)a = Bx{k)ax + By{k)ay + 
Bz{k)az (an overall factor oi g^B has been absorbed into 
the definition of B to simplify notation). The directions 
of Bint as dependent on the directions of k are shown 
in FigP. The SIA-type term gives a Bint that keeps 
its orientation as k{\\E) is rotated between the [110] and 
[110] directions, while both BIA-type Bint coming from 
Hi and i?3 change their sign between [110] and [110]. 
The difference between Hi and H^ is that the latter has 
a finite Bint when E\\x,y whereas the former has zero 
Bint for the same directions. 
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[110] 
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H3 : BIA Strain Split 



FIG. 1: Direction of the internal magnetic field Bi„t for the 
three spin splitting Hamiltonians Hi, H2, H3, {\,Ci,D > 0) 
considering the electric field E (and hence the average mo- 
mentum) to be in plane. It was also assumed that txx = (-yy, 
as appropriate in the experiment The experimental SIA 
data cannot be explained by the term H2, but nevertheless, 
the SIA-type Bint seen in the experiment has the same di- 
rection as the one plotted here. When -E||a; or E\\y, the BIA 
term induces a Bint\\E whereas the BIA term H\ does not 
induce spin splitting for these directions. This constitutes a 
simple check of the experiment 



In the values of the splitting (3 are measured on 
the [110] and [lIO] directions and because of the sign- 
changing properties of Bint mentioned above, the BIA 
and the SIA contributions to (3 can be obtained as follows: 

f^BiA = (/3[iio] - /3[iio])/2, PsiA = (/3[iio] +/3[iio])/2. 
Surprisingly, the spin splitting is more of a BIA-type 
rather than an SIA-type, contradicting the conventional 
knowledge that an SIA-type term described by H2 is re- 
sponsible for s pin splitting in strained semiconductors 

|2lllll2llilii&. 

Theoretically, the Dresselhaus term Hi is a bulk- 
inversion asymmetry term that appears even in the ab- 
sence of strain. As observed in the experiment, the ficti- 
tious internal magnetic field Bint is perpendicular to the 
momentum k: Bint(k)k = 0, where Bx = Xkx{ky — kf), 
By and Bz being obtained by cubic permutation. For 
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GaAs, the constant A = 22eVA^. However, we beheve 
this term is not responsible for the spin sphtting observed 
in the experiment The observed sphtting is hnear in 
momentum fc, inconsistent with the Hi ~ k^. Exper- 
iments performed on InSb |29l| . another material with 
inversion asymmetry, support this conclusion and point 
strongly to the fact that strained InSb is described by 
H2. In [i^ stress of up to 4 kbar is applied mechan- 
ically on a 1 X 1 X lOmnT' sample and Shubnikov-de- 
Haas oscilations are used to probe the band structure. 
Without applied strain, the conduction band exhibits 
a spin-splitting that is small and cubic in fc, described 
by Hi. In 2^ the application of diagonal does not in- 
duce any observable spin splitting whereas the applica- 
tion of shear strain induces a splitting linear in fc, de- 
scribed by H2. Relatively large stress- induced splitting 
of the Fermi surfaces occurs in the lower concentration 
{n = 1.4 X 10^5 - 2.0 X 10"cm-3) samples H^. The 
energy splitting dispersion switches from in the un- 
strained case to k when strain (stress) above Ikbar is ap- 
plied, in accordance to H2 becoming dominant over Hi. 
From a theoretical estimate, at n = 3 x lO^^cm"'^, Hi 
should be of the same size as H3 and roughly one order 
of magnitude lower than H2. The spin splitting at the 
Fermi wavevector kp = 0.96 x 10^m~^ due to Hi is less 
than 10~^eV. By contrast the spin splitting due to H2 is 
C^exykp = 5.06 x IQ-'^t^ycV = 5 x lO^^ - 15 x IQ-'^eV 
for txy = 0.1% — 0.3% as in [l| (see Table[l] for conversion 
of strain components from [l| to the orthogonal system 
txy). An experimental value of Ca/fi. = 8 x lO^m/s for 
GaAs was used "SO^. Contrary to previous remarks fsij, 
there is hence no theoretical or experimental a-priori rea- 
son to disregard the strain-dependent spin splitting terms 
in favor of the Dresselhaus k^ term for the doping values 
in the experiment 1]. 

There are a number of experimental reasons inLJ hint- 
ing the marginal significance of the k^ term. In [l| strain 
plays a critical role in generating the spin-orbit coupling 
Bint. Samples prepared from the same wafer but un- 
strained show a reduction by an order of magnitude of 
Bint along both the [110] and the [110] directions. If Hi 
were responsible for spin splitting, its value would remain 
unchanged upon varying strain. Strain could only enter 
the system through the variation of the effective electron 
mass in the x, y, z directions, as ;31] points out. However, 
these variations with strain are of a maximum 2% — 3% 
[23. thereby not accounting for the order of magni- 
tude variation of the spin-splitting between the strained 
and the unstrained cases observed in [J]. 

The term H2 is a structural inversion asymmetry 
(SIA)-type term that has its origin in the acoustic phonon 
interaction of the valence band with the conduction band 
|23j ]. In the framework of the Kane's 8x8 matrix (2x2 
for the conduction and split-off band and 4x4 for the 
valence band) the conduction band couples to the va- 
lence band. In systems with inversion symmetry where 



the selection rules for L are satisfied, it is impossible 
to couple spin-0 (|s)) with spin-1 (\z)) through a spin- 
2 term (e.y) and hence {s\txy\z) — 0. However, when 
inversion symmetry is broken, the fore-mentioned term 
need not be zero as the L selection rule need not apply. 
Upon straining, the matrix elements between the con- 
duction and valence band have the form (sjea^yjz) (plus 
cyclic permutations) where \s) is the s-orbital and j^) 
is one of the p orbitals. Through perturbation theory, 
one can compute the effect of this valence-conduction 
band interaction when projected to the conduction band 
and obtain the conduction band effective Hamiltonian 
H2 HHMmiii]. Taking into account that the electric 
field is in-plane (< k^ >= 0) and that Cxy 7^ (see 
Table[l]), in H2 the components of the internal mag- 
netic field (which due to the rescaling by giiB has units 
of energy) are: B^ = ^Csexyky, By = -^C^eyxK _= 
— ^Cz^xykx- Switching coordinates to the [110] and [110] 
directions, -B[iio] — ^C^e^yk^^nQ]^, -Sj^jg] = —^C^exyk^uo] 
(see Figpp. Since H2 is an SIA term, the spin splitting 
(3 will be of SIA type /J^'}^ (th stands for the theoretical 
estimate). Since < k >= \Tnvd where Vd is the drift ve- 
locity of the spin packed due to the electric field, we get 
a simple formula for the 

PfiA = ^^.y^ (2) 

By using the experimentally known value for C^/h — 
8 X lO^m/2, the predicted values for /Sg/^ are given in 
Table |3 The theoretical values are larger than the ob- 
served ones by a factor of 3 — 30 and no matching trend 
between the data and the SIA term H2 can be found. 
Moreover, as remarked in Q no systematic correlation 
between the experimentally observed SIA contribution 
and the strain is observed. We hence come to the con- 
clusion that the SIA spin splitting observed in is not 
induced by the uniform shear strain (which would give 
the values /3gj^ and which have been confirmed in me- 
chanical experiments) but borrows substantially from the 
dislocations and strain gradient inherent in growing such 
a thick sample through MBE techniques. This is not to 
say that the SIA term is negligible: as seen in Table |3 
the SIA term is substantial and comparable in magni- 
tude with the BIA term. However, the SIA term does 
not correlate with strain ad cannot be described by H2. 

The remaining spin splitting term is H^. Although 
this term is allowed by group theory, it only shows up 
at higher order in perturbation theory than H2 in the 
k ■ p method. We claim that in the experiment Q] this 
term is responsible for the spin-splitting observed, and 
determine the value of the constant D. We note, how- 
ever, that this does not settle the theoretical puzzle of 
why the H3 would be more significant than H2 in this 
case, which might have to do with the conditions of 
the experiment such as low temperature, the appear- 
ance of dislocations and strain gradient. H3 is a Hamil- 
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Sample 




^xx — ^yy 




riexp 
PSIA 


l3flA 


r,exp 

Pbia 




D/h 


A 


0.46 


-0.16 


0.2 


-24 


604 


75 


121 


1.59 


B 


0.14 


-0.2 


0.08 


-26 


241 


13 


38 


0.5 


E 


0.13 


-0.42 


0.18 


69 


543 


43 


78 


1.03 


F 


0.07 


-0.32 


0.12 


54 


362 


31 


79 


1.04 


G 


0.04 


-0.32 


-0.04 


44 


-121 


31 


86 


1.13 


H 


0.13 


-0.42 


0.26 


-2 


785 


24 


43 


0.56 


I 


0.04 


-0.16 


0.2 


65 


604 


23 


115 


1.51 



TABLE I: Strain components (in %) in 1] converted to 
cartesian coordinates (C and D are compositionally graded 
films Q). The conversion equations, obtained after a sim- 
ple coordinate transformation, are: e^z ~ e[ooi] , ^xx ~ ^yy ~ 
i(<;[iio] + £[110]), = |(<;[iioi - £[110])- The measured spin 
splitting slope values Pg^/j^ and as well as the theoret- 

ical value /Sgj^ obtained using the Hamiltonians H2 for the 
SIA - type splitting are given in neV nS jj.m~^ , the units of 
Q. The value of the deformation potential D in the term 
H3, determined as explained in Eq(|lJ is given for each sam- 
ple. The discrepancies between the values of and /3sj^ 
make the case that H2 is not the full mechanism for strain 
SIA splitting. For five (A, E, F, G, I) out of the seven sam- 
ples the values of the deformation potential D deduced from 
the H3 model of the strain BIA splitting are within 30% of 
each other. Combined with the simplicity of the model, this 
match gives weight to the suggestion that BIA splitting is 
produced by a BIA term. ~ ^^x) is measured in 
lO^neV" nS fim ^, while D is measured in 10'*m/s 

tonian of BIA-type and vanishes with vanishing strain, 
thereby satisfying two of the experimental observations 
in P]. From Eq.pQ, the internal magnetic field Bint reads: 

-E^x — Dkx{f^zz ^yy)' -^y — ^^yi^xx ^xx^: Bz — 

Dkzicyy — txx), and since the in-plane electric field in- 
fluences only the in-plane momentum, < >= and 
only the two in-plane components of the internal mag- 
netic field remain. In accordance with the experiment, 
we place e^x — ^yy, and e^z — £xx > (Table[l]). We 
hence have B.^ = k.xD{e^^ - e^x), By = kyD{exx - <^zz), 
D{ezz—ixx) > 0. Since i3i„t is not perpendicular to fc| 
Bxkx + Byky = {kl - kl)D{tzz - <^xx), one may think 
this term is incompatible with the observed Bint -!l E\\k 
in This, however, would be hasty: the experiment 
is performed in only two directions, with ^^[[[llO] and 
i?||[110], for which kx = ^ky. For these two directions 
only, the Bint in Hj, is perpendicular to the momentum 
and the electric field, hence satisfying a major constraint 
the experimental data poses on the theory. Since the 
value of the constant D is unknown from previous exper- 
imental studies (although it was suggested that they can 
be sometimes sizable |33 ) there is no way of theoretically 
predicting the values of the spin-splitting from our model. 
However, we can check if the model is consistent with the 
experimental data and we can also obtain a value of the 
constant D which, being a material constant, should be 
similar on all the samples cited here. Since < k >= y-Vd 



where Vd is the spin drift velocity along the spin packet 
we find: 

f3[BIA]^^ji^zz-^xx)m (3) 

We can determine the value of D from the experimental 
data for /3 and strain e: 

D ^ 1 P[BIA] 

h 2m tzz - exx 

As a consistency check, since 13 is a material constant, 
P[B^iA\/{.^zz — ^xx) should be quasi-constant between the 
samples quoted in the experiment. In 5 out of the 7 
samples studied in [ij, the values of P[B'iA]/i^^^ ~ ^^^) 
are close together to within 30%, lumped in two groups 
(samples A, I are very close to each other, and within 
30% of the value for E, F, G which are again very close 
between themselves). The samples E, F, G were grown in 
the same day. The deviant samples B, H were also grown 
in the same day, and hence the variation of the constant 
coefficient D within a sample set that was grown on the 
same day is less than 15% (21|. Different growth con- 
ditions are most likely responsible for the (still small) 
variations between samples grown in different days. The 
consistency check is further proof that H3 is the term re- 
sponsible to the spin-splitting in . The values obtained 
for D are given in Table p. 

We showed that H3 is a BIA-type Hamiltonian vanish- 
ing with vanishing strain, with an internal magnetic field 
that is perpendicular to the applied electric field for the 
two experimental directions [110] and [110] and which is 
consistent with the reported data for the spin splitting. 
On the other hand. Hi and H2, the previously known 
spin splitting terms, fail to reproduce the data on more 
than several counts. It is easy to experimentally prove, 
using the setup in , that H3 is responsible for the spin 
splitting is easy: one would measure the internal mag- 
netic field due to BIA on the x or y direction. In this 
case, an H3 term would give an internal magnetic field 
parallel to E (of course, there will also be an internal B 
from an SIA term that is still perpendicular to E, but 
a component of Bint parrallel to the electric field should 
be easily detectable). 

In another beautiful experiment, Kato et al. measure 
through Farraday Rotation (FR) a nonzero uniform mag- 
netization Pel induced by driving an electric current (elec- 
tric field) through the sample E of their previous experi- 
ment . It has been long predicted 1^ that semicon- 
ductors with spin-orbit coupling will exhibit a uniform 
magnetization when placed in an electric field generating 
a charge current. This can be trivially understood by a 
simple argument: writing the spin-orbit Hamiltonian as 
a k dependent magnetic field Zeeman coupled to spin, 
Bint{k)(T{k), the application of an electric field E will 
make the average value of the momentum be non-zero 
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< k >= ^Et where r is the momentum relaxation time. 
This creates a non-zero average < Bint >= B{< k >) 
which orients the spins along its direction through the 
Zeeman-like coupling of the spin-orbit term. 

We now try to numerically estimate the value of the 
uniform magnetization p^i using the BIA-type H3. From 
0, the BIA contribution to the uniform magnetization 

can be obtained as p^i^^ = ^(pf — ^"^^ 
around 3 x lO^^m^^ for E ~ W^V/m. We will now try to 
estimate this from first principles using Hg, as the main 
BIA term and using the value of D for sample E deduced 
in Table[l]. A simple linear response calculation of the 
magnetization (Ti to the electric current Jj (due to the 
applied electric field Ej) gives: 




27rer 



Pel 



Qij = {TcJiJj) = 



h 

d^k nE_ 



■QijEj 



{2nf 



dB OB, 
B,— -B- 



dki 



dki 



(5) 



where i,j = x,y, z and Bi{k) are the components of the 
internal magnetic field for iJs, B — 



n E± is the Fermi function of the spin-split energies E± — 
■^k^ ± B for H3. For the Hamiltonian H3, considering 
frrr — ^yy wc obtaiu for: 



Pel 



BIA 



eTmkFD{tz 



-E 



(6) 



As previously pointed out the magnetization is paral- 
lel to the electric field for E\\x or E\\y. This provides 
an important and easy check of the above assumption 
that the observed strain spin splitting comes from H3. 
The only two directions where Pei is perpendicular to the 
electric field are [110] and [110], the directions on which 
the experiment is performed. Considering a sample of 
mobility fj, = 0.6m^/Vs 

we obtain an estimate for 
p^i — 3.45 X lO^^m"^ for a field E — lO^V/m, compared 
to an experimental value of 3 x 10^*to~'^ for the same 
value of the electric field. The theoretical value obtained 
is within the experiment's error margins. 

Finally, using the current setup in [ij, ^ we propose an 
experiment to test the prediction of dissipationless spin 
current. For spin 1/2 two-dimensional systems, the ini- 
tial prediction is subject to some sort of controversy, 
|14L 133 . 134I] as the introduction of impurities apparently 
makes the spin current vanish. We here adopt the al- 
ternative view and propose a clear-cut experiment which 
can see the spin accumulation due to the spin current. 
Similar to the 2D case, in the present case, the appli- 
cation of an electric field Ej to a semiconductor with 
spin orbit coupling will create a spin current flowing 
perpendicular to the electric field and polarized perpen- 
dicular to both the field and the direction of flow. Using 



FIG. 2: An electric field will cause both an observed net uni- 
form bulk magnetization and a proposed spin current. The 
spin current will accumulate at the edges over a spin diffu- 
sion length of more than 500nm thus making its detection 
practical with a beam slightly more focused than in 2]. The 
uniform magnetization and the spin current spin accumulatio 
are 7r/2 out of phase. 



linear response, the expression for the spin conductance 
is: 



Ji - <^ijEj 
fi2 f (fk nE_ - nE+ dB„, 



2m J (27r)3 



where i, j,l,m,n = x,y,z, tinm is the totally antisym- 
metric tensor in 3 dimensions and Bi{k) are the com- 
ponents of the internal magnetic held. For TJs and for 
^xx = ^yy the only non-zero components of the spin con- 
ductance are: 



1 



f (fik nE_ - riE^ 2 



^21 



'12 



hD{e,,-eyy)mJ (2^)3 (fc^ fc2)3/2 



1 



2m (27r) 



27r 



d(j) / desin^ (j){k^_-kl) 



(8) 



where cf), 9 are the polar angles of k and where k- , k^ are 
the fermi momenta of the two bands. When both bands 
are occupied (positive Fermi energy), we find — fc+ = 
X{k) = D{ezz - exx)ksm6. Usually the spin 
splitting is much smaller than the Fermi energy, and we 
can define an average Fermi momentum fcj? ~ ^i^- + 
k+) « (37r^n)^/3^ being the dopant density. With this, 
we find that the spin-conductivity will be independent of 
the value of the strain: 



'21 



'12 



(9) 
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The result for the spin conductance is intermediate be- 
tween the 2D spin 1/2 spin current and the 3D spin 3/2 
spin current. Similar to but unlike |^ the spin conduc- 
tance depends on the fermi momenturn, a characteristic 
of the 3D. Unlike 0, but similar to the spin con- 
ductance does not depend on the strength of spin-orbit 
coupling. Even though the spin conductance does not de- 
pend on the value of strain, it is essential that spin-orbit 
splitting (due to strain in this case) be present. Upon the 
application of an electric field on the x axis, a spin current 
will flow on the y axis spin polarized in the z direction. 



For n 
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16, 



and a field E = IQ'^V/ra we esti- 
mate a spin current jspin = Jij:^^ = 2 x 10^^(/XB/cm^s) 
where is a Bohr magneton. Since spin conductivity 
varies as na and charge conductivity varies as n, for low 
values of n the spin conductance will overtake the charge 
conductance and the spin current will be larger than the 
charge current caused by the electric field. The density 

-.2/3 / 2e (67r^l^ 



where /i is 
for a 



at which this happens is n ''^ < 
the mobility in the sample, or n < 2 x 
sample of mobility ^ — Q.Qrn^ /Vs. 

The flow will result in accumulation on the opposite zx 
faces of the crystal (see Fig [2]). For the present experi- 
ment, we estimate this spin accumulation of the order of 
spin Ts = lO^^^B /cm^. Due to the extremely spin life 
time of above Ins, the distance from the edge of the sam- 
ple, the spin diffusion length is very large, of the order 
L = 500nm — 1/im. The FR beam used in j2| has a resolu- 
tion of 4.7/im/9.7/im on the x and y axis respectively, but 
focusing the beam within l/im is possible [iSIl^- Then, 
if the spin current prediction is right, applying the FR 
beam on the edge of the sample should give a clear sig- 
nal (larger than the uniform magnetization in the bulk). 
Since the uniform magnetization and the spin accumula- 
tion due to spin hall current are perpendicular to each 
other, in time-resolved FR experiments, the spin current 
spin accumulation and the uniform magnetization are out 
of phase by 7r/2 (see Figj^J 

In conclusion, we have analyzed two very recent exper- 
iments P, and proved that conventional spin splitting 
terms and strain spin splitting terms do not explain the 
data. We have introduced a previously largely unknown 
term and made the case as to why it explains the observed 
features in 0,01 . We have proposed further simple exper- 
iments to verify our assertions. If true, our proposal gives 
rise to the clear possibility of obtaining a uniform magne- 
tization parallel to the applied electric field, as opposed 
to the one perpendicular to it that has been observed so 
far. Along with predicting a 3D spin current, we have 
also proposed a way to test the spin currents in spin 1/2 
systems. 
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